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ABSTRACT 

We discuss the effects of an isolated gravitational lens on the rotation of the plane of polar- 
ization of linearly polarized light rays, the so called gravitational Faraday rotation, in metric 
theories of gravity. By applying the thin lens approximation, we derive simple expressions 
for extended, rotating astrophysical systems deflecting electromagnetic radiation from back- 
ground sources. Higher order corrections and the case of a Kerr black hole are also considered 
to show how the rotation of polarization is a purely gravito-magnetic effect. Prospects of a fu- 
ture detection of a gravitational Faraday rotation are examined. 
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1 INTRODUCTION 

Conceivable post-Newtonian theories of gravity, such as general 
relativity, predict space-time curvature by mass-energy currents rel- 
ative to other masses. This feature is known as intrinsic gravito- 
magnetism l Ciufolini & Wheeler 1995 l. 

The occurrence of several phenomena is due to gravito- 
magnetism. Small test-gyroscopes rotate with respect to distant 
stars in the vicinity of a spinning body due to its rotation. This 
is called the Lense-Thirring precession. Together with test parti- 
cles and gyroscopes, also photons and clocks probe some gravito- 
magnetic phenomena. Bending and time delay of electromagnetic 
waves are two important phenomena predicted by theories of grav- 
ity. Gravito-magnetic effects due to the spin of an astrophysical ob- 
ject might be detected in gravitational lensing events l Sereno 2002 
ISereno & C ardone 2002 1. The rotation of the lens induces particu- 
lar signatures on both the position and magnification of multiple 
images of a background source, also perturbing the microlensing 
induced amplification curve I Sereno 2003b l. 

Another well-known gravitational effect due to space-time 
curvature by mass currents is the rotation of the plane of polar- 
ization for linearly polarized light rays. Under geometric optics, 
a light ray follows a null geodesic regardless of its polarization 
state and the polarization vector is parallel transported along the 
ray ( Misner 'et al. 19731 . Such a gravitational rotation of the plane 
of polarization in stationary space-times is a gravitational ana- 
logue of the electro-magnetic Faraday effect iPiran & Safier 1985 
Ishihara et al. 1988, Nouri-Zonoz 1999 1, i.e., the rotation that a 
light ray undergoes when passing through plasma in the presence of 
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a magnetic field. The gravitational effect is known as gravitational 
Faraday rotation or Rytov or Skrotskii effect. 

From the original investigations in Skrotskii 119571 and 
Rytov 119381 . the effect of a gravitational field on the plane of 
polarization of electromagnetic waves has been addressed by 
several authors tBalazs 19581 IPlebanski 19rjQl IGodfrey 1970 
ISu & Mallett 19801 IConnors et al. 19801 IPiran & Safier 1985 
llshihara et al. 19881 |Dyer & Shaver T992) INouri-ZonozT999 
Kopeikin & Mashoon 2002 , Sereno 2004 1. The case of a weak 
gravitational field has deserved particular attention. In Pleban- 
ski 119601 . the Maxwell's equations in the gravitational field 
of an isolated physical system were solved and the rotation of 
the plane of polarization around the propagation vector was 
considered. In Kopeikin & Mashoon 1 2002 1, a formula describing 
the Skrotskii effect for arbitrary translational and rotational motion 
of gravitating bodies was derived by solving the equations of 
motion of a light ray in the first post-Minkowskian approximation. 
In Sereno 1 2004 1, the general formula for the angle of rotation of 
the plane of polarization of a linearly polarized electromagnetic 
wave in a stationary space-time was re-obtained with an heuristic 
approach based on Mach's principle on dragging of inertial frames, 
without integrating the equations of motion. 

The assumption of the existence of a dynamical space-time 
curvature, as opposed to a flat space-time of special relativity, is 
shared among different viable theories of gravity, such as gen- 
eral relativity, Brans-Dicke theory and the Rosen bimetric the- 
ory. Metric theories of gravity are defined such that l Wil l 19931 
Ciufolini & Wheeler 1995 l: i) space- time is a Lorentzian manifold; 
ii) the world lines of test bodies are geodesies; Hi) the equivalence 
principle holds in the medium strong form. The usual rules for the 
motion of particles and photons in a given metric still apply, but the 
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metric may be different from that derived from the Einstein's field 
equations. 

A comparison among general relativity and other viable the- 
ories of gravity can be led with suitable tests on the basis of 
higher-order effects, such as intrinsic gravito-magnetism. To date, 
results from laser-ranged satellites provide the only demonstration 
of the gravito-magnetic field. In 1995-2002, the Lense-Thirring 
precession, due to the Earth spin, was measured by studying 
the orbital perturbations of LAGEOS and LAGEOS II satellites 
(Ciufol ini & Pavlis 1998 1. Its experimental value agrees with about 
20% accuracy with the prediction of general relativity. The NASA's 
Gravity Probe-B satellite should improve this measurement to an 
accuracy of 1%. 

In this paper, we discuss the gravitational Faraday ro- 
tation by an isolated rotating gravitational deflector in the 
standard framework of gravitational lensing l Pette rset al. 200 fl 
ISchneider et al. 19 92 1. when the source of radiation and the ob- 
server are remote from the source of the gravitational field (lens). 
As known,thanks to some usual approximations, the effect of 
the spin of the deflector on the lensing potential can be treated 
in a quite simple way ( Sereno 2002 1. The lens equation can be 
modified to consider how some observable quantities, like posi- 
tion and amplification of images, are perturbed by the rotation 
( Sereno & Cardone 2002 , Sereno 2003b l. Here, we want to extend 
such a formalism to the Skrotskii effect and to consider some ex- 
tended lenses of astrophysical interest. Standard assumptions allow 
us to compute higher-order approximation terms in the calculation, 
so that a general treatment of the Skrotskii effect can be performed 
in metric theories of gravity. 

The paper is as follows. In Sec.|2| the weak-field, slow motion 
approximation in metric theories of gravity is introduced and the 
weak field limit of the gravitational Faraday rotation is performed. 
In Section [3] we derive some simple formulae, using the thin lens 
approximation, to treat systems of astrophysical interest. Higher or- 
der correction to the Skrotskii effect are discussed in Sec.|4j then, 
the case of light rays propagating in the vacuum region outside a 
Kerr black hole is considered. Section [5] contains some final con- 
siderations. 



2 THE WEAK FIELD LIMIT 



Standard hypotheses of gravitational lensing iPette rset al. 200 ll 
ISchneider et al. 1992i assume that the gravitational lens is local- 
ized in a very small region of the sky and its lensing effect is weak. 
In Sereno ( 2003a l, an approximate metric element generated by an 
isolated mass distribution was written in the weak field regime and 
slow motion approximation, up to the post-post-Newtonian order, 
and with non-diagonal components accounting for effects of grav- 
ity by currents of mass. In an asymptotically Cartesian coordinate 
system, the metric g a p can be expressed as 



1 + 2^- + 0( £ 4 )) c 2 dt 2 - (l - 2 7 | + G( 



quantifies the contribution to the space-time curvature of the mass- 
energy currents. It measures the strength of the intrinsic gravito- 
magnetic field l Ciufolini & Wheeler 1995 l. In general relativity, 
/i = 1; in Newtonian theory, fi = 0. Additional terms in a 
parametrized expansion of the metric element can also be con- 
sidered (M atzner et al. 19811 IRichter & Matzner 1982i . Preferred 
frame effects, violations of conservation of four momentum, pre- 
ferred location effects are not considered in our approximate metric 
element. In Eq. Q, (j> is the Newtonian potential, 
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where G is the gravitational Newton constant and p is the mass 
density of the deflector; <?i/c 2 is order ~ 0(e 2 ). 
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V I =V + V +0{e ) ) is a vector potential taking into ac- 
count the gravito-magnetic field produced by mass currents. To the 
lowest order of approximation, 



Vi (t,x)~-G 
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where v is the velocity field of the mass elements of the deflector 2 . 

We assume that, during the time light rays take to traverse the 
lens, the potentials in Eqs. ^2|[3} vary negligibly little. Then, the 
lens can be treated as stationary. In Eqs. ^2||3j. we have neglected 
the retardation I Schneider et al. 1992 1. 

The polarization vector is dragged along by the rotation of the 
inertial frames. Once we have introduced the notation 



h = goo, Ai = , 

3oo 
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the net angle of rotation around the unit tangent three-vector k 
along the path between the source and the observer, in a station- 
ary space- time, reads iNouri-Zonoz 1999 Sereno 2004 1 



fiRk = — — 



VhVx Adx, 



(5) 



where dx = kdlp and dip is the spatial distance in terms of the 
spatial metric jij (Landau & Lifshits 1985 1, 
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gij + ^iSM j dx^daP = jijdx t dx :i 
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where e <C 1 denotes the order of approximation. In the 
above metric element, 7 is a standard coefficient of the 
post-Newtonian parametrized expansion of the metric tensor 
( Ciufo lini & Wheeler 1995IIWill 1993L measuring space curvature 
produced by mass. In general relativity, it is 7 = 1; in the Brans- 
Dicke theory, 7 = j^j. M is a non standard parameter which 



1 Latin indices run from 1 to 3, whereas Greek indices run from to 3. 



Operations on three-vectors are defined in the three-dimensional 
space with metric 7^. 

In Dyer & Shiver 1 1992 1, it was shown, using symmetry ar- 
guments in a relativistic context, that, for astrophysically interest- 
ing cases of both observer and source at large distances from the 
lens, a nonshperical, nonrotating lens cannot cause a net rotation of 
the polarization vector. This argument is confirmed by Eq. : the 
Skrotskii effect is due to mass-currents, i.e., it is a purely gravito- 
magnetic phenomenon. 

. ^E. (fyHx^J \SdWf r \tY' m ^1" me orc ^ er °f approximation is de- 
termined by off-diagonal components of the metric. The main term 
is order ~ C(e 3 ). To determine a gravito-magnetic effect to the 
second leading order, i.e., to 0(e 5 ), we need to consider terms in 
goo and ga up to 0(e 2 ). Higher order terms would give contribu- 
tions ~ C(e 4 ) x V~ 0{e 7 ) to the Faraday rotation, flsk is at best 

2 Bold symbols denote spatial vectors. The controvariant components of 
spatial three- vectors are equal to the spatial components of the correspond- 
ing four- vectors. 
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~ C(e 3 ) unlike the bending of light, which is ~ 0{e 2 ). That is 
why the gravitational Faraday rotation is usually neglected in light 
propagation analyses. 

In the weak field limit, h and A are simply related to the grav- 
itational potentials. It is 



A, 



1 + 2^ + 0(e% 



the proper arc length reads 

dh ~ |l- 7 l+C( e 3 )|d/ E , 



(7) 
(8) 

(9) 



imation assumes that rays of electromagnetic radiation propagate 
along straight lines. The integration along the line of sight (l.o.s.) is 
accurate enough to evaluate the main contribution to the Skrotskii 
effect. To this order, we will employ the unperturbed Minkowski 
metric r/ a p = (1, — 1, — 1, —1) for operations on vectors. 

Let us define the weighted average of a quantity Q along the 
line of sight, 



where dlv = y Sijdx x dxi is the Euclidean arc length, whereas the 
spatial metric reduces to 

7 ^ = (l-2^+0( e 4 ))<%, (10) 

To calculate the gravitational Faraday rotation to order G N , 
we need the path of the deflected light ray to the order G N ~ 1 
I Fishbach & Freeman 1980 Sereno 2003a l. Taking the a; 3 -axis 
along the unperturbed photon path, to the first order in G, the unit 
propagation vector reads k(i) ~ {C(e 2 ), C(e 2 ), 1 + 70, }• The 
Faraday rotation to order 0(e 5 ) turns out 
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where comma stands for differentiation and p is the spatial projec- 
tion of the null geodesies. The first term in the right hand side of 
Eq. QTJ 

gives the main contribution to the gravitational Faraday 
rotation l Sereno 2004 1. Other terms represent higher-order correc- 
tions. The second and the third contributions account for the inte- 
gration along the deflected path p. The fourth term in Eq. <l It de- 
rives from the difference between the gravito-magnetic potentials 
A and V. 



3 THE THIN LENS APPROXIMATION 

In almost all astrophysical systems, the gravitational lens can be 
treated as geometrically thin. The extent of the lens in the direction 
of the incoming ray is small compared to both the distances be- 
tween lens and observer and lens and source, so that the maximal 
deviation of the actual ray from a light path propagating through 
unperturbed space-time is small compared to the length scale on 
which the gravitational field changes. In this picture, gravitational 
effects are quite local to the lens vicinity. 

The lens plane is transverse to the incoming, unperturbed light 
ray direction. The position vector in the lens plane is the spatial vec- 
tor £ = {x 1 , x 2 , O} . We introduce the coordinate I along the x 3 - 
axis, so that the lens plane corresponds to I — 0. The Bom approx- 
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where E is the surface mass density of the deflector, 
E«)= / p(U)dl. 



(12) 
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The leading contribution to the Skrotskii effect can be written 
; <Sereno 20041 
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Inserting the expression for the gravito-magnetic potential and per- 
forming the integration along the line of sight, we get 
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where Li. .$. is the angular momentum, per unit mass, of a mass 
element at £ with respect a photon with impact parameter £, 

(£i.o. s .)i.o.s. = (A. .s.(£ 2 - £ 2 ') - (A-o.s.^ 1 - e 1 '); (i6) 

v 1 and v 2 are the components of v along the and the £ 2 -axes, 
respectively. 

The angular momentum of the lens can determine a direction 
of rotation for the polarization vector. When its projection along the 
line of sight is not null, a handness is introduced in the system and 
some rotation can be caused. Once the deflector is oriented such 
that (I/i. .s.)i.o.s. = 0, the photon moves like in an equatorial plane 
and a mirror symmetry about the lens plane holds. Such a symmetry 
cancels out the effect (Dyer & Shaver 1992) . 

3.1 Spherical lenses 

In most of the astrophysical situations, lenses can be approximated 
with nearly spherical mass density profiles. Let us consider spher- 
ical deflector in rigid rotation (constant angular velocity w), when 
v = tiJXx. We limit to a slow rotation so that the deformation 
caused by rotation is negligible and the body has a spherical sym- 
metry. 

Taking the centre of the source as the spatial origin of a back- 
ground inertial frame, the weighted average of the rotational veloc- 
ity along the line of sight turns out to be 

(v)i.o.s. = {C 2 ^ 1 _ C 1 ^ 2 ; — £, 2 ui.o. s ., £ Wi.o.s.} • (17) 

We can, now, evaluate the integral in Eq. <15t ; the polarization vec- 
tor of a linearly polarized light ray with impact parameter |£| is 
rotated by 



n. 
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where Mc y i(> |£|) is the projected lens mass outside |£| 
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The gravitational Faraday rotation is proportional to the mass out- 
side |£| and can be significant for lenses with slowly decreasing 
mass density. f2sk can be compared to the contribution to bending 
of light due to the gravito-magnetic field, SIgrm (see Eqs. (15, 16) 
in Sereno & Cardone ( 2002 1). Typically, Q,s\a and Qgrm are of the 
same order. Contrary to the Faraday effect, the gravito-magnetic 
contribution to the deflection angle depends only on the component 
of the angular momentum in the lens plane ( Sereno 2002 1. Whereas 
the components of the angular momentum along the £ l -axis, related 
to the momentum of inertia of the mass within £ about a central 
axis, contribute only to Hgrm I Sereno & Cardone 2002 1, the con- 
tribution due to the mass outside the impact parameter appears in 
both flsk and JIgrm- 

Let us consider an isothermal sphere (IS). The IS is a density 
profile widely used to model systems on very different scales, from 
galaxy haloes to clusters of galaxies. The surface mass density is 



1 



2G yiiF+^F' 



(20) 



where <t v is the velocity dispersion and £ c is a finite core radius. 
Since the total mass is divergent, we introduce a cut-off radius R 3> 
|£|. The gravitational Faraday rotation is 



) = ^aW a . s . (y&T&- vWTg) 



(21) 



In particular, in the inner regions (|£| <C R), the above equation 
reduces to 



^m«R)^^aWo. s R. 



(22) 



In the case of lensing of distant quasars by foreground galax- 
ies, images may form inside the galaxy radius. We can model a 
typical galaxy as a singular (£ c = 0) IS with a v ~ 200 km s , 
R ~ 20 Kpc and J ~ 0.1 Mq Kpc s , as derived from nu- 
merical simulations I Vitviska et al. 2002 1. From Eq. I22i . we find 

— 6 milliarcsec. 

Let us consider Ogrm f° r an isothermal sphere. It is 
<Sereno & Cardone 20021 . 



For a source at rest at the origin of the coordinates, assuming 
that the polar axis of the coordinate system coincides with the rota- 
tion axis, the weak field limit of the Kerr metric in Boyer-Lindquist 
coordinates reads 



ds ~ dsa 



4:llG J . 2/1/ , w \ 

— 3 a sin o(rd<p)(cdt), 



(24) 



where ds| ch is the line element for a spherically symmetric solu- 
tion and where terms of quadratic and higher order in the angular 
momentum J have been neglected. 

The metric in Eq. 1241 can be expressed in an equivalent 
isotropic form by introducing a new radius variable, p, such as 



GM 



(25) 



where M is the total mass of the black hole. Substituting Eq. <25t 
in Eq. 1241 and introducing quasi-Minkovskian coordinates related 
to (p, 8, <fi) by the usual transformation rules, we get the weak field 
limit of the Kerr metric in the isotropic form. By arbitrarily orien- 
tating the angular momentum J, we get 



ds 2 
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where x = \J 5ijX l x : > and 
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In the above equations and in what follows, operations on vectors 
are performed using rj a /3. Eqs. I27II28I are also valid for large dis- 
tances from an isolated physical system. Since, 



„ s GJ-3(J-x)x 
Vx V (x) = L— 
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Slgk and r^GRM are °f tne same order but depend on different com- 
ponents of the angular momentum of the lens. 



4 HIGHER ORDER CORRECTIONS 

To determine higher order terms in the gravitational Faraday rota- 
tion, it is enough to consider only the monopole term in the gravito- 
electric field but we must go over the dipole moment in the gravito- 
magnetic potential. We want now discuss higher order terms in a 
rotating black hole. 



4.1 The Kerr black hole 

The Faraday effect in a Kerr metric, in the weak field limit, has 
been already discussed by several authors llshihar aet al. 19881 
INouri- Zonoz 1999 1. The right order of approximation in the cal- 
culation has been ascertained, but a disagreement on the numerical 
value still persists. Now, we want to apply our formalism to light 
rays passing through the vacuum region outside a rotating black 
hole. 



it is easy to verify that, to order C(e 3 ), there is no gravitational 
Faraday rotation l Pleba nski 19601 |Kopeikin & Mashoo n" 2002| 
|Godfrey 1970| . Let us consider higher order terms. 
The unit wave-vector, to order C(e 2 ) reads 



k ( i) 




(30) 



where iisch = 2GM/c 2 is the Schwarzschild radius of the lens. 
The deflected path to order G is actually the path in a spherically 
symmetric metric. To this order, the photon path lies in the plane 
through source, lens and observer. 

After lengthy but straightforward calculations, we can perform 
the integrations in Eq. <1 H . We find 3 



Kerr T G M Jl.o.s. 7\ 

fisk ~-p-— —+0(e) 



(31) 



3 Our Ji.o.a. corresponds to —am cos #o in Ishihara et al. 1198 8 1 and 
Nouri-Zonoz 1 1999 i. 
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As already known l Ishihara et al. 1988 , Nouri-Zonoz 1999 1, the net 
rotation is of order C(e 5 ). 

Only the projection of the angular momentum along the line- 
of-sight enters the effect. In Nouri-Zonoz 1 1999 1, the integration 
is performed along the unperturbed path. However, contributions 
from the deflected path cancel out. The parameter 7, which mea- 
sures space curvature produced by mass, does not enter the effect. 



5 CONCLUSIONS 

We have discussed the theory of the gravitational Faraday rotation 
in the weak field limit in a metric theory of gravity. Due to a gravito- 
magnetic field, a net rotation of the polarization plane can occur. 
Unlike the phenomena of bending of light and time delay of elec- 
tromagnetic waves, gravito-magnetism appears in the leading term 
of the rotation angle of the plane of polarization. 

Useful formulae have been derived for relevant astrophys- 
ical systems. By applying the thin lens approximation, we ob- 
tained quite compact expressions for extended spinning gravita- 
tional lenses. The gravitational Faraday increases for lenses with 
a slowly decreasing mass density. The effect of a gravito-magnetic 
field is of the same order in both bending of light and rotation of the 
plane of polarization. Contrary to the gravito-magnetic deflection, 
the Skrotskii effect depends only on the component of the angu- 
lar momentum along the line of sight. Furthermore, whereas the 
rotation of the deflector can break the cylindrical symmetry of the 
deflection angle in the lens plane, no azimuthal dependence enters 
the gravitational Faraday rotation. A non null projection of the total 
angular momentum of the lens along the line of sight introduces a 
handness in the system but no preferred direction in the alignment 
of the light beam. 

To consider higher order effects, we have considered light rays 
propagating outside a rotating black hole. Due to the order of ap- 
proximation, the effect on light rays propagating in the vacuum re- 
gion outside the event horizon of a Kerr black hole has been often 
missed. Applying a method based on the use of the Walker-Penrose 
constant up to higher-order terms, Ishihara et al. 1 1988 1 first calcu- 
lated the angle of rotation due to the presence of the black hole's 
spin. Nouri-Zonoz 1 19991 used the 1+3 formulation of stationary 
space-times and the quasi-Maxwell form of the vacuum Einstein 
equations to study Kerr spaces. The results in Ishihara et al. 1 1988 1 
and Nouri-Zonoz 1 1999 1 are of the same order but differ for a nu- 
merical factor. Here, we have properly considered all the contribu- 
tions to the gravitational Faraday rotation, In particular, the effect 
of the curved path has been accounted for. The gravitational Fara- 
day rotation depends linearly on fj,, the parameter which quantify 
and measures the gravito-magnetic effect. Even at higher order of 
approximation, ji is the only parameter which enters the effects. 
Other parametrized post-Newtonian parameters, such as 7, do not 
appear. For fj, = 1, our result agrees with Nouri-Zonoz 1 1999 1. 

Techniques involving the angular relationship between the 
intrinsic polarization vectors and the morphological structure 
of extended radio jets have been proposed to probe the 
mass distribution of foreground lensing galaxies iKronbe rg 1991 1 
Kronberg et al. 1996 1. The method has been then extended in Surpi 

6 Harari 119991 to consider the effect of weak-lensing by large 
scale structure. The basic principle underlying these investigations 
is that the polarization vector of the radiation would not be rotated 
by the gravitational field of the deflector. 

Here, we have considered gravitational Faraday rotation by as- 
trophysically significant lens models, such as an isothermal sphere. 



The Skrotskii effect turns out to be negligible with respect to today 
observational uncertainties. Anyway, high quality data in total flux 
density, percentage polarization and polarization position angle at 
radio frequencies already exist for multiple images of some gravi- 
tational lensing systems, like B02 18+357 (Biggs et al. 1999 1. Since 
the impressive development of technological capabilities, it is not a 
far hypothesis to obtain observational evidences in a next future. 

Furthermore, polarization measurements of the cosmic mi- 
crowave background radiation should probe very subtle effects, 
such as primordial gravity wave signals or a remnant from the 
reionization epoch. A full understanding of the whole foreground 
contamination turns out to be crucial in data analysis. So, a proper 
modeling and quantification of the gravitational Faraday effect is 
necessary. 
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